3.2.8 \(\int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx\) [108]

3.2.8.1 Optimal result
3.2.8.2 Mathematica [B] (verified)
3.2.8.3 Rubi [A] (verified)
3.2.8.4 Maple [B] (verified)
3.2.8.5 Fricas [F(-1)]
3.2.8.6 Sympy [F]
3.2.8.7 Maxima [F]
3.2.8.8 Giac [F]
3.2.8.9 Mupad [F(-1)]

3.2.8.1 Optimal result

Integrand size = 37, antiderivative size = 228 \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\frac {2 \sqrt {-d g+c h} (a+b x) \sqrt {\frac {(b g-a h) (c+d x)}{(d g-c h) (a+b x)}} \sqrt {\frac {(b g-a h) (e+f x)}{(f g-e h) (a+b x)}} \operatorname {EllipticPi}\left (-\frac {b (d g-c h)}{(b c-a d) h},\arcsin \left (\frac {\sqrt {b c-a d} \sqrt {g+h x}}{\sqrt {-d g+c h} \sqrt {a+b x}}\right ),\frac {(b e-a f) (d g-c h)}{(b c-a d) (f g-e h)}\right )}{\sqrt {b c-a d} h \sqrt {c+d x} \sqrt {e+f x}} \]

output
2*(b*x+a)*EllipticPi((-a*d+b*c)^(1/2)*(h*x+g)^(1/2)/(c*h-d*g)^(1/2)/(b*x+a 
)^(1/2),-b*(-c*h+d*g)/(-a*d+b*c)/h,((-a*f+b*e)*(-c*h+d*g)/(-a*d+b*c)/(-e*h 
+f*g))^(1/2))*(c*h-d*g)^(1/2)*((-a*h+b*g)*(d*x+c)/(-c*h+d*g)/(b*x+a))^(1/2 
)*((-a*h+b*g)*(f*x+e)/(-e*h+f*g)/(b*x+a))^(1/2)/h/(-a*d+b*c)^(1/2)/(d*x+c) 
^(1/2)/(f*x+e)^(1/2)
 
3.2.8.2 Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(583\) vs. \(2(228)=456\).

Time = 31.71 (sec) , antiderivative size = 583, normalized size of antiderivative = 2.56 \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=-\frac {2 \sqrt {\frac {(d g-c h) (a+b x)}{(b g-a h) (c+d x)}} (c+d x)^{3/2} \left (\frac {a d \sqrt {\frac {(d g-c h) (e+f x)}{(f g-e h) (c+d x)}} (g+h x) \operatorname {EllipticF}\left (\arcsin \left (\sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}\right ),\frac {(b c-a d) (-f g+e h)}{(d e-c f) (b g-a h)}\right )}{(d g-c h) (c+d x) \sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}}+\frac {b c \sqrt {\frac {(d g-c h) (e+f x)}{(f g-e h) (c+d x)}} (g+h x) \operatorname {EllipticF}\left (\arcsin \left (\sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}\right ),\frac {(b c-a d) (-f g+e h)}{(d e-c f) (b g-a h)}\right )}{(-d g+c h) (c+d x) \sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}}+\frac {b (f g-e h) \sqrt {\frac {(-d e+c f) (d g-c h) (e+f x) (g+h x)}{(f g-e h)^2 (c+d x)^2}} \operatorname {EllipticPi}\left (\frac {d (-f g+e h)}{(d e-c f) h},\arcsin \left (\sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}\right ),\frac {(b c-a d) (-f g+e h)}{(d e-c f) (b g-a h)}\right )}{(d e-c f) h}\right )}{d \sqrt {a+b x} \sqrt {e+f x} \sqrt {g+h x}} \]

input
Integrate[Sqrt[a + b*x]/(Sqrt[c + d*x]*Sqrt[e + f*x]*Sqrt[g + h*x]),x]
 
output
(-2*Sqrt[((d*g - c*h)*(a + b*x))/((b*g - a*h)*(c + d*x))]*(c + d*x)^(3/2)* 
((a*d*Sqrt[((d*g - c*h)*(e + f*x))/((f*g - e*h)*(c + d*x))]*(g + h*x)*Elli 
pticF[ArcSin[Sqrt[((-(d*e) + c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]], (( 
b*c - a*d)*(-(f*g) + e*h))/((d*e - c*f)*(b*g - a*h))])/((d*g - c*h)*(c + d 
*x)*Sqrt[((-(d*e) + c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]) + (b*c*Sqrt[ 
((d*g - c*h)*(e + f*x))/((f*g - e*h)*(c + d*x))]*(g + h*x)*EllipticF[ArcSi 
n[Sqrt[((-(d*e) + c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]], ((b*c - a*d)* 
(-(f*g) + e*h))/((d*e - c*f)*(b*g - a*h))])/((-(d*g) + c*h)*(c + d*x)*Sqrt 
[((-(d*e) + c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]) + (b*(f*g - e*h)*Sqr 
t[((-(d*e) + c*f)*(d*g - c*h)*(e + f*x)*(g + h*x))/((f*g - e*h)^2*(c + d*x 
)^2)]*EllipticPi[(d*(-(f*g) + e*h))/((d*e - c*f)*h), ArcSin[Sqrt[((-(d*e) 
+ c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]], ((b*c - a*d)*(-(f*g) + e*h))/ 
((d*e - c*f)*(b*g - a*h))])/((d*e - c*f)*h)))/(d*Sqrt[a + b*x]*Sqrt[e + f* 
x]*Sqrt[g + h*x])
 
3.2.8.3 Rubi [A] (verified)

Time = 0.32 (sec) , antiderivative size = 228, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.054, Rules used = {183, 412}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx\)

\(\Big \downarrow \) 183

\(\displaystyle \frac {2 (a+b x) \sqrt {\frac {(c+d x) (b g-a h)}{(a+b x) (d g-c h)}} \sqrt {\frac {(e+f x) (b g-a h)}{(a+b x) (f g-e h)}} \int \frac {1}{\left (h-\frac {b (g+h x)}{a+b x}\right ) \sqrt {\frac {(b c-a d) (g+h x)}{(d g-c h) (a+b x)}+1} \sqrt {\frac {(b e-a f) (g+h x)}{(f g-e h) (a+b x)}+1}}d\frac {\sqrt {g+h x}}{\sqrt {a+b x}}}{\sqrt {c+d x} \sqrt {e+f x}}\)

\(\Big \downarrow \) 412

\(\displaystyle \frac {2 (a+b x) \sqrt {c h-d g} \sqrt {\frac {(c+d x) (b g-a h)}{(a+b x) (d g-c h)}} \sqrt {\frac {(e+f x) (b g-a h)}{(a+b x) (f g-e h)}} \operatorname {EllipticPi}\left (-\frac {b (d g-c h)}{(b c-a d) h},\arcsin \left (\frac {\sqrt {b c-a d} \sqrt {g+h x}}{\sqrt {c h-d g} \sqrt {a+b x}}\right ),\frac {(b e-a f) (d g-c h)}{(b c-a d) (f g-e h)}\right )}{h \sqrt {c+d x} \sqrt {e+f x} \sqrt {b c-a d}}\)

input
Int[Sqrt[a + b*x]/(Sqrt[c + d*x]*Sqrt[e + f*x]*Sqrt[g + h*x]),x]
 
output
(2*Sqrt[-(d*g) + c*h]*(a + b*x)*Sqrt[((b*g - a*h)*(c + d*x))/((d*g - c*h)* 
(a + b*x))]*Sqrt[((b*g - a*h)*(e + f*x))/((f*g - e*h)*(a + b*x))]*Elliptic 
Pi[-((b*(d*g - c*h))/((b*c - a*d)*h)), ArcSin[(Sqrt[b*c - a*d]*Sqrt[g + h* 
x])/(Sqrt[-(d*g) + c*h]*Sqrt[a + b*x])], ((b*e - a*f)*(d*g - c*h))/((b*c - 
 a*d)*(f*g - e*h))])/(Sqrt[b*c - a*d]*h*Sqrt[c + d*x]*Sqrt[e + f*x])
 

3.2.8.3.1 Defintions of rubi rules used

rule 183
Int[Sqrt[(a_.) + (b_.)*(x_)]/(Sqrt[(c_.) + (d_.)*(x_)]*Sqrt[(e_.) + (f_.)*( 
x_)]*Sqrt[(g_.) + (h_.)*(x_)]), x_] :> Simp[2*(a + b*x)*Sqrt[(b*g - a*h)*(( 
c + d*x)/((d*g - c*h)*(a + b*x)))]*(Sqrt[(b*g - a*h)*((e + f*x)/((f*g - e*h 
)*(a + b*x)))]/(Sqrt[c + d*x]*Sqrt[e + f*x]))   Subst[Int[1/((h - b*x^2)*Sq 
rt[1 + (b*c - a*d)*(x^2/(d*g - c*h))]*Sqrt[1 + (b*e - a*f)*(x^2/(f*g - e*h) 
)]), x], x, Sqrt[g + h*x]/Sqrt[a + b*x]], x] /; FreeQ[{a, b, c, d, e, f, g, 
 h}, x]
 

rule 412
Int[1/(((a_) + (b_.)*(x_)^2)*Sqrt[(c_) + (d_.)*(x_)^2]*Sqrt[(e_) + (f_.)*(x 
_)^2]), x_Symbol] :> Simp[(1/(a*Sqrt[c]*Sqrt[e]*Rt[-d/c, 2]))*EllipticPi[b* 
(c/(a*d)), ArcSin[Rt[-d/c, 2]*x], c*(f/(d*e))], x] /; FreeQ[{a, b, c, d, e, 
 f}, x] &&  !GtQ[d/c, 0] && GtQ[c, 0] && GtQ[e, 0] &&  !( !GtQ[f/e, 0] && S 
implerSqrtQ[-f/e, -d/c])
 
3.2.8.4 Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(847\) vs. \(2(209)=418\).

Time = 1.29 (sec) , antiderivative size = 848, normalized size of antiderivative = 3.72

method result size
elliptic \(\frac {\sqrt {\left (b x +a \right ) \left (d x +c \right ) \left (f x +e \right ) \left (h x +g \right )}\, \left (\frac {2 a \left (\frac {g}{h}-\frac {a}{b}\right ) \sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \left (x +\frac {c}{d}\right )^{2} \sqrt {\frac {\left (-\frac {c}{d}+\frac {a}{b}\right ) \left (x +\frac {e}{f}\right )}{\left (-\frac {e}{f}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \sqrt {\frac {\left (-\frac {c}{d}+\frac {a}{b}\right ) \left (x +\frac {g}{h}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, F\left (\sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}, \sqrt {\frac {\left (\frac {e}{f}-\frac {c}{d}\right ) \left (\frac {g}{h}-\frac {a}{b}\right )}{\left (-\frac {a}{b}+\frac {e}{f}\right ) \left (-\frac {c}{d}+\frac {g}{h}\right )}}\right )}{\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (-\frac {c}{d}+\frac {a}{b}\right ) \sqrt {b d f h \left (x +\frac {a}{b}\right ) \left (x +\frac {c}{d}\right ) \left (x +\frac {e}{f}\right ) \left (x +\frac {g}{h}\right )}}+\frac {2 b \left (\frac {g}{h}-\frac {a}{b}\right ) \sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \left (x +\frac {c}{d}\right )^{2} \sqrt {\frac {\left (-\frac {c}{d}+\frac {a}{b}\right ) \left (x +\frac {e}{f}\right )}{\left (-\frac {e}{f}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \sqrt {\frac {\left (-\frac {c}{d}+\frac {a}{b}\right ) \left (x +\frac {g}{h}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \left (-\frac {c F\left (\sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}, \sqrt {\frac {\left (\frac {e}{f}-\frac {c}{d}\right ) \left (\frac {g}{h}-\frac {a}{b}\right )}{\left (-\frac {a}{b}+\frac {e}{f}\right ) \left (-\frac {c}{d}+\frac {g}{h}\right )}}\right )}{d}+\left (\frac {c}{d}-\frac {a}{b}\right ) \Pi \left (\sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}, \frac {-\frac {g}{h}+\frac {a}{b}}{-\frac {g}{h}+\frac {c}{d}}, \sqrt {\frac {\left (\frac {e}{f}-\frac {c}{d}\right ) \left (\frac {g}{h}-\frac {a}{b}\right )}{\left (-\frac {a}{b}+\frac {e}{f}\right ) \left (-\frac {c}{d}+\frac {g}{h}\right )}}\right )\right )}{\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (-\frac {c}{d}+\frac {a}{b}\right ) \sqrt {b d f h \left (x +\frac {a}{b}\right ) \left (x +\frac {c}{d}\right ) \left (x +\frac {e}{f}\right ) \left (x +\frac {g}{h}\right )}}\right )}{\sqrt {b x +a}\, \sqrt {d x +c}\, \sqrt {f x +e}\, \sqrt {h x +g}}\) \(848\)
default \(\text {Expression too large to display}\) \(1250\)

input
int((b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2),x,method=_RETU 
RNVERBOSE)
 
output
((b*x+a)*(d*x+c)*(f*x+e)*(h*x+g))^(1/2)/(b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e 
)^(1/2)/(h*x+g)^(1/2)*(2*a*(g/h-a/b)*((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/(x+c/d 
))^(1/2)*(x+c/d)^2*((-c/d+a/b)*(x+e/f)/(-e/f+a/b)/(x+c/d))^(1/2)*((-c/d+a/ 
b)*(x+g/h)/(-g/h+a/b)/(x+c/d))^(1/2)/(-g/h+c/d)/(-c/d+a/b)/(b*d*f*h*(x+a/b 
)*(x+c/d)*(x+e/f)*(x+g/h))^(1/2)*EllipticF(((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/ 
(x+c/d))^(1/2),((e/f-c/d)*(g/h-a/b)/(-a/b+e/f)/(-c/d+g/h))^(1/2))+2*b*(g/h 
-a/b)*((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/(x+c/d))^(1/2)*(x+c/d)^2*((-c/d+a/b)* 
(x+e/f)/(-e/f+a/b)/(x+c/d))^(1/2)*((-c/d+a/b)*(x+g/h)/(-g/h+a/b)/(x+c/d))^ 
(1/2)/(-g/h+c/d)/(-c/d+a/b)/(b*d*f*h*(x+a/b)*(x+c/d)*(x+e/f)*(x+g/h))^(1/2 
)*(-c/d*EllipticF(((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/(x+c/d))^(1/2),((e/f-c/d) 
*(g/h-a/b)/(-a/b+e/f)/(-c/d+g/h))^(1/2))+(c/d-a/b)*EllipticPi(((-g/h+c/d)* 
(x+a/b)/(-g/h+a/b)/(x+c/d))^(1/2),(-g/h+a/b)/(-g/h+c/d),((e/f-c/d)*(g/h-a/ 
b)/(-a/b+e/f)/(-c/d+g/h))^(1/2))))
 
3.2.8.5 Fricas [F(-1)]

Timed out. \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\text {Timed out} \]

input
integrate((b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2),x, algor 
ithm="fricas")
 
output
Timed out
 
3.2.8.6 Sympy [F]

\[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\int \frac {\sqrt {a + b x}}{\sqrt {c + d x} \sqrt {e + f x} \sqrt {g + h x}}\, dx \]

input
integrate((b*x+a)**(1/2)/(d*x+c)**(1/2)/(f*x+e)**(1/2)/(h*x+g)**(1/2),x)
 
output
Integral(sqrt(a + b*x)/(sqrt(c + d*x)*sqrt(e + f*x)*sqrt(g + h*x)), x)
 
3.2.8.7 Maxima [F]

\[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\int { \frac {\sqrt {b x + a}}{\sqrt {d x + c} \sqrt {f x + e} \sqrt {h x + g}} \,d x } \]

input
integrate((b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2),x, algor 
ithm="maxima")
 
output
integrate(sqrt(b*x + a)/(sqrt(d*x + c)*sqrt(f*x + e)*sqrt(h*x + g)), x)
 
3.2.8.8 Giac [F]

\[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\int { \frac {\sqrt {b x + a}}{\sqrt {d x + c} \sqrt {f x + e} \sqrt {h x + g}} \,d x } \]

input
integrate((b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2),x, algor 
ithm="giac")
 
output
integrate(sqrt(b*x + a)/(sqrt(d*x + c)*sqrt(f*x + e)*sqrt(h*x + g)), x)
 
3.2.8.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\int \frac {\sqrt {a+b\,x}}{\sqrt {e+f\,x}\,\sqrt {g+h\,x}\,\sqrt {c+d\,x}} \,d x \]

input
int((a + b*x)^(1/2)/((e + f*x)^(1/2)*(g + h*x)^(1/2)*(c + d*x)^(1/2)),x)
 
output
int((a + b*x)^(1/2)/((e + f*x)^(1/2)*(g + h*x)^(1/2)*(c + d*x)^(1/2)), x)